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This  study  focuses  on  optimum  operating  strategies  for  liquid-fed  direct  methanol  fuel  cells  (DMFCs)  to 
minimize  methanol  consumption.  A  mathematical  model  is  developed  and  verified  with  experimental 
data  from  the  literature  using  the  parameter  estimation  method.  The  model  consists  of  a  set  of  differen¬ 
tial  and  algebraic  equations  and  makes  it  possible  to  describe  zero  initial  hold-up  conditions.  Based  on  the 
model,  steady-state  simulation  results  are  obtained  and  explain  the  dependence  on  the  feed  concentration 
of  key  variables  such  as  cell  voltage,  cell  power  density,  overpotentials  of  both  electrodes,  and  methanol 
crossover  ratio.  Dynamic  simulation  results  are  also  presented  to  check  the  transient  behaviour  of  a  DMFC 
operated  from  start-up  to  shut-down.  Dynamic  optimization  allows  determination  of  the  optimum  tran¬ 
sient  strategies  of  feed  concentration  required  to  maximize  the  fuel  efficiency.  With  six  scenarios  of  power 
density  load,  it  is  demonstrated  that  the  optimum  transient  strategies  depend  heavily  on  both  the  load  of 
power  density  and  the  number  of  control  actions.  The  main  advantage  of  these  approaches  is  to  reduce 
fuel  consumption  and,  ultimately,  to  enable  DMFCs  to  be  operated  more  efficiently. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  is  a  hot  research  topic 
because  it  represents  a  power  source  for  portable  electronic  and 
vehicular  applications,  employs  a  solid  polymer  electrolyte  that 
works  at  relatively  low  temperature,  and  uses  an  inexpensive  liquid 
fuel  with  an  especially  high  specific  energy.  Methanol  has  a  the¬ 
oretical  energy  density  of  6.1  Whg-1.  Even  a  20%  efficient  DMFC 
is  more  than  three  times  efficient  than  advanced  lithium  batter¬ 
ies  [1].  Nevertheless,  in  practice,  DMFC  suffers  from  two  intrinsic 
drawbacks.  One  is  the  slow  reaction  kinetics  of  the  methanol 
oxidation  in  the  catalyst  layer  of  the  negative  electrode  (anode) 
because  of  the  complexity  of  its  mechanism.  This  results  in  a 
much  lower  open-circuit  voltage  (OCV)  than  is  seen  with  other 
types  of  fuel  cells,  such  as  the  hydrogen  proton-exchange  mem¬ 
brane  fuel  cell  (PEMFC),  and,  consequently,  a  far  lower  power 
generation  for  a  given  size.  The  other  drawback  is  that  methanol 
permeates  into  the  proton-exchange  membrane  (PEM)  to  reach 
the  cathode  catalyst  layer,  where  it  forms  a  mixed  potential  that 
decreases  the  cell  voltage.  This  phenomenon,  termed  methanol 
crossover,  not  only  results  in  a  waste  of  fuel  but  also  degrades 
the  performance  of  the  fuel  cell  [2].  It  is  thus  very  important 
to  predict  methanol  crossover  and  to  find  operating  strategies 
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that  can  reduce  methanol  crossover,  so  as  to  improve  system  effi¬ 
ciency. 

Recently,  extensive  studies  have  been  carried  out  to  predict 
methanol  crossover  in  DMFCs;  such  studies  have  employed  not  only 
experimental  methods  [1-7]  but  also  modelling  systems  [8-15]. 
Danilov  et  al.  [8]  developed  a  three-dimensional,  two-phase  CFD 
model  to  predict  the  internal  phenomena  of  a  DMFC.  Sundmacher 
et  al.  [9]  and  Zhou  et  al.  [10]  formulated  a  dynamic  differential  and 
algebraic  equations  (DAEs)-based  model  and  simulated  cell  volt¬ 
age  response  to  dynamic  changes  of  methanol  feed  concentration 
in  order  to  determine  whether  periodically  pulsed  methanol  feed¬ 
ing  could  reduce  methanol  crossover.  These  authors  used  an  anodic 
reaction  mechanism  that  consisted  of  four  steps,  and  derived  an 
analytical  solution  of  anodic  reaction  rate,  assuming  the  first  reac¬ 
tion  step  to  be  rate-determining.  Jeng  and  Chen  [11]  developed  a 
one-dimensional  steady-state  model  on  the  anode  side  of  the  DMFC 
and  found  that  methanol  crossover  became  serious  at  low  current 
density  and  high  methanol  concentration.  Schultz  and  Sundmacher 
[12]  improved  on  their  previous  work  [9,10]  with  a  partially  one¬ 
dimensional  dynamic  model  to  characterize  the  flow  directional 
behaviour  in  the  PEM  and  diffusion  layers.  Krewera  and  Sund¬ 
macher  [13]  applied  transfer  function  analysis  to  investigations  of 
the  dynamic  behaviour  of  DMFCs.  They  investigated  cell  voltage 
responses  to  cell  current  step  changes,  and  developed  a  criterion 
for  cell  voltage  overshooting. 

Another  approach  to  modelling  is  dynamic  optimization  to  find 
the  optimum  operating  conditions,  which  was  first  addressed  by  Xu 
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Nomenclature 

As  cross-sectional  electrode  area  (m2) 

q  concentration  of  component  i  in  the  anode  compart¬ 

ment  (molm-3) 

cf*  concentration  of  component  i  in  the  anode  catalyst 
layer  (molm-3) 

cfeed  feed  concentration  of  component  i  (mol  m-3 ) 

Cj  double  layer  capacity  of  the  electrode  j  ( F  m-2 ) 

dM  thickness  of  the  membrane  (m) 

Dj^JeOH  diffusion  coefficient  of  methanol  in  the  membrane 

(m2  s-1) 

F  Faraday’s  constant  (96,485  C  mol-1 ) 

zceii  cell  current  density  (A  m-2 ) 

kj  reaction  rate  constant  of  the  electrode  j 

(molm-2  s-1) 

kLS  mass  transfer  coefficient  (m  s-1 ) 
n?M  mass  flux  of  component  i  in  the  membrane 

(molm-2  s-1) 

Pj  pressure  of  the  electrode  j  (Pa) 

rj  reaction  rate  of  the  electrode  j  (mol  m-2  s-1 ) 

R  universal  gas  constant  (8.314J  mol-1  K-1 ) 

t  time  (s) 

T  temperature  (I<) 

Uce\\  cell  voltage  (V) 

U standard  cell  voltage  (1.21  V) 
vM  superficial  velocity  of  water  in  the  membrane 
(ms-1) 

Va  volume  of  the  anode  compartment  (m3 ) 

V F  anode  feed  flow  rate  (m3  s-1 ) 

V£l  volume  of  the  anode  catalyst  layer  (m3 ) 

xfL  mole  fraction  of  component  i  in  the  anode  catalyst 
layer 

Greek  letters 

oij  charge  transfer  coefficient  of  the  electrode  j 

£ void  fraction  of  catalyst  layer  of  the  electrode  j 
Pj  overpotential  of  the  electrode  j  (V) 

/cM  conductivity  of  the  membrane  (£2-1  m-1 ) 

Subscripts 
a  anode 

c  cathode 

C02  carbon  dioxide 

MeOH  methanol 

Superscripts 
CL  catalyst  layer 

Feed  feed 

M  membrane  (PEM) 

STD  at  standard  conditions 


et  al.  [14].  They  maximized  the  integration  of  cell  voltage  over  the 
entire  time  domain  under  constant  cell  current  density  and  adopted 
the  model  of  Sundmacher  et  al.  [9]  to  obtain  the  optimum  feed  con¬ 
centration.  More  practical  optimization  work  was  performed  by  Ko 
et  al.  [15],  who  developed  a  non-isothermal  dynamic  optimization 
model  based  on  the  isothermal  model  of  Sundmacher  et  al.  [9]  and 
optimized  the  fuel  efficiency  satisfying  the  power  load.  According  to 
their  results,  of  the  three  operating  variables,  namely,  the  methanol 
concentration,  feed  flow  rate  at  the  anode,  and  air  flow  rate  at  the 
cathode,  the  feed  concentration  has  the  strongest  influence  on  the 
dynamic  behaviour  of  a  DMFC. 


In  this  research,  we  focus  on  the  identification  of  the  optimum 
operating  strategies  of  DMFCs  over  their  entire  operating  proce¬ 
dure,  including  start-up  operation.  In  the  initial  state  of  the  start-up 
operation,  the  cell  is  assumed  to  be  fully  hydrated;  it  contains  nei¬ 
ther  methanol  nor  carbon  dioxide.  In  order  to  simulate  and  optimize 
the  start-up  operation,  it  is  necessary  to  consider  zero  initial  hold¬ 
up  conditions  for  methanol  and  carbon  dioxide  in  the  dynamic 
model  of  DMFC,  even  though  these  conditions  tend  to  cause  unex¬ 
pected  numerical  problems.  In  contrast  to  previous  contributions 
to  the  dynamic  optimization  of  DMFCs  [14,15],  the  present  study 
adopts  zero  initial  hold-up  conditions  to  describe  the  early  state  of 
the  start-up  operation. 


2.  Mathematical  model 

The  focus  of  the  research  is  a  single  cell  of  a  DMFC.  As  shown  in 
Fig.  1,  the  cell  consists  of  a  membrane  electrode  assembly  (MEA) 
and  anode  and  cathode  flow-fields  with  a  symmetric  structure  that 
centres  on  a  PEM.  The  MEA  can  be  divided  into  five  layers:  a  PEM, 
two  catalyst  layers,  and  two  diffusion  layers.  The  PEM  separates  the 
anode  side  from  the  cathode  side;  ideally,  it  allows  proton  transport 
only.  A  dilute  solution  of  methanol  (typically  less  than  2  M  for  the 
active  type,  and  about  5  M  for  the  passive  type)  is  fed  to  the  anode 
flow  channel;  air  or  pure  oxygen  is  supplied  to  the  cathode  side. 
The  feeds  permeate  the  diffusion  layers  and  reach  the  catalyst  lay¬ 
ers.  At  the  anode  catalyst  layer,  1  mol  of  methanol  reacts  with  an 
equimolar  amount  of  water  to  produce  1  mol  of  carbon  dioxide  and 
6  mol  of  protons  and  electrons.  By  contrast,  at  the  cathode  catalyst 
layer,  1.5  mol  of  oxygen  react  with  6  mol  of  protons  and  electrons 
to  produce  3  mol  of  water. 


2.1.  Assumptions  and  conditions 

The  aim  is  to  develop  a  rigorous  dynamic  model  to  describe  and 
analyze  the  dynamic  behaviour  of  a  DMFC.  The  model  is  a  modi¬ 
fication  of  those  of  Sundmacher  et  al.  [9]  and  Zhou  et  al.  [10],  and 
considers  material  transports  and  electrochemical  reactions  under 
the  following  assumptions  and  conditions. 


•  The  system  is  treated  as  a  continuous  stirred  tank  reactor  (CSTR) 
that  is  operated  isothermally. 

•  The  diffusion  layers  and  the  flow  fields  are  combined  into  the 
compartments. 

•  Ohmic  drops  in  current-collectors  and  electric  connections  are 
negligible. 

•  Oxygen  is  supplied  in  excess  and  therefore  no  oxygen  mass  bal¬ 
ance  is  required. 

•  The  PEM  is  fully  hydrated  and  allows  for  transport  of  protons  and 
methanol  only. 

•  The  water  concentration  is  constant  (excess  component  in  liquid 
mixture). 

•  Mass  transport  resistances  in  the  catalyst  layers  are  negligible 
because  these  are  thin  in  comparison  with  the  diffusion  layers 
and  PEM. 

•  The  mass  transport  coefficients  of  methanol  and  carbon  dioxide 
are  equal. 

•  The  electrochemical  reactions  of  the  anode  and  the  cathode  cat¬ 
alyst  layers  follow  the  Butler-Volmer  reaction  mechanism. 

•  The  mixtures  in  the  system  are  treated  as  pure  liquid  phase 
mixtures,  i.e.,  gas-phase  formation  by  release  of  carbon  dioxide 
bubbles  is  not  considered. 
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]  Compartment  (diffusion  layer  +  flow  channel) 

|  Catalyst  layer  (anode:  carbon-supported  Pt-Ru,  cathode:  carbon-supported  Pt) 
!  V>  Proton  exchange  membrane  (Nation®  117) 


Anode:  CHflH{l)  +  H20(l)  PtRu  >  C02(g)  +  6//+  +  6e~  =  0.02  V 

Cathode:  |02(g)  +  6//+  +  6e~  — >  3 H20(l)  =123  V 

Overall:  CH .OH (/)  +  j02(g)  - >  C02(g)  +  2H 20(l)  U*t„  =  1.21  K 

Fig.  1.  Basic  structure  and  reaction  scheme  of  DMFC. 


2.2.  Model  equations 


anode  and  cathode  sides  are  therefore  given  by 


In  the  anode  compartment,  methanol  is  fed  to  reach  the  anode 
catalyst  layer,  and  carbon  dioxide  comes  from  the  anode  catalyst 
layer.  The  material  balances  of  methanol  and  carbon  dioxide  in  the 
anode  compartment  can  be  written  as 


=  T(ice„  -  /a)  =  T(iCell  -  6 Fra)  (7) 

^  =  T(-icell  -  ic)  =  T(_jcel|  _  6 Frc  +  4Fn“e0H)  (8) 


dCMeOH 

df 

dcco2 

dt 


“  y  ^CMeOH  cMeOH) 


kisAs 


(CMeOH  CMeOH 


^F/rFeed  r  \  klSAS  ,  CL  * 

Y~  (  cco2  _  Cc°2 )  w  ^  Cc°2  ~  CC°2  ' 


Assuming  that  the  PEM  is  fully  hydrated  and  the  total  flux 
through  the  PEM  is  almost  the  same  as  the  flux  of  water  (only 
because  the  concentration  of  methanol  is  extremely  low  in  the 
DMFC),  Jeng  and  Cheng  [11  ]  derived  the  rate  of  methanol  crossover 
as  follows: 


In  the  anode  catalyst  layer,  methanol  from  the  anode  com¬ 
partment  is  consumed  by  methanol  crossover  as  well  as  by 
electrochemical  reaction  producing  equimolar  carbon  dioxide: 


dcCL 

ULMeOH 

df 


/<lsAs 

~W 


As 

(QvieOH  -  CMeOF[)  -  ^CL^MeOH  +  ra) 


dc§32  kisAs  CL  ^S 

dt  ~~  yCL  ^C02  CC02)+yCLra 


(3) 

(4) 


It  is  well  known  that  the  mechanism  of  the  electrochemical  reac¬ 
tions  in  a  DMFC  is  complicated  and  utilizes  several  intermediate 
reactions  [9,10].  To  reduce  the  number  of  unknown  parameters, 
however,  this  research  adopts  Butler-Volmer-type  equations  for 
the  rate  expressions  of  methanol  oxidation  at  the  anode  catalyst 
layer  and  oxygen  reduction  at  the  cathode  catalyst  layer,  i.e., 


ra  =  ka  (x^oh  eXP  ( ' ~  XC02  eXP  ("  ^  ^  ^))  (5) 

rc  =  kc  ((^)15  exp("Sr,,c)  ■exp((ir,c))  (6) 

where  xfL  is  the  mole  fraction  of  component  i  (methanol,  car¬ 
bon  dioxide  and  water)  and  p02  is  the  partial  pressure  of  oxygen 
in  the  cathode  catalytic  layer.  At  the  anode  side,  six  protons 
and  electrons  are  produced  from  the  methanol  oxidation;  at 
the  cathode  side,  four  protons  and  electrons  are  produced  from 
the  undesired  methanol  oxidation  due  to  methanol  crossover 
(CH3OH  +  O.5O2  C02  +4H+  +4e-),  whereas  six  protons  and  elec¬ 

trons  are  consumed  by  oxygen  reduction.  The  charge  balances  of 


„M  =  4LeOHVMexP(vMdM/D^eOH)  m 

Me0H  (exp(vMdM/Dj^e0H)  -  1 ) 

where  vM  is  the  superficial  velocity  of  water  in  the  PEM,  which  can 
be  calculated  by 


vM  =  MH2Q  ^H2Q  (  /cell  A 

Ph2o  \  F  ) 


(10) 


The  overall  cell  voltage  is  a  function  of  the  open-circuit  voltage 
(OCV),  anode  and  cathode  overpotentials,  and  ohmic  loss  in  the 
PEM: 


Uceii  =  u™  -riz  +  ric-  ^i'ceii  (H ) 

When  the  cell  is  operated  under  galvanostatic  conditions,  the  over¬ 
all  cell  voltage  can  be  calculated  if  only  two  overpotentials  are 
obtained. 

To  simulate  the  cell  using  Eqs.  ( 1 )-( 11 ),  it  is  necessary  to  find  sev¬ 
eral  physical  properties  especially  on  the  PEM  (Nation®  117).  They 
are  usually  expressed  as  a  function  of  temperature.  The  diffusivity 
of  methanol  in  the  PEM,  Dj^[e0H  in  Eq.  (9),  is  known  to  be  a  function 
of  temperature  [16],  namely: 

D^eOH  =  ^eOH.ref  ?XP  (  jT  -  }))  (12) 

Scott  et  al.  [16]  reported  that  Springer  et  al.  [17]  obtained  a  value 
of  2436  K  for  A E/R  that  pertains  to  measurements  on  Nation®  117. 
For  a  reference  value,  Kauranen  and  Skou  [18]  measured  a  superfi¬ 
cial  diffusivity  of  methanol  in  the  PEM  of  4.9  x  10-10  m2  s-1  at  60  °C. 
Finally,  the  diffusivity  of  methanol  in  the  PEM  at  temperature  T  can 
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Excess  Air  H20(g) 


be  calculated  by: 

D“e0H  =4.9  X  10-10  exp  (2436  -  1))  (13) 

The  temperature  dependence  of  the  electro-osmotic  drag  coef¬ 
ficient  of  water  (Ah2o  in  Eq.  (10))  in  fully  hydrated  PEM  can  be 
estimated  from  the  experimental  data  of  Ren  et  al.  [  19]  by  regression 
analysis  using  an  exponential  growth  function  with  three  parame¬ 
ters: 

Xh2o  =  -0.0382  +  0.1585  exp(8.7129  x  1CT3T)  (14) 

Springer  et  al.  [17]  obtained  the  ionic  conductivity  of  Nation® 
117  as  a  function  of  temperature  and  water  content  ijr. 

km  =  (0.5139t/f  -  0.326)  exp  (l268  (^  -  1))  (15) 

The  water  content  (ifr)  of  the  membrane  is  defined  as  the  ratio  of 
the  number  of  water  molecules  to  the  number  of  S03~  groups  per 


Table  1 

Operating  conditions  and  parameters 


Parameter 

Value 

Operating  conditions 

Anode  inlet  pressure,  Pa  (Pa) 

101,325 

Cathode  inlet  pressure,  Pc  (Pa) 

303,975 

Operating  Temperature,  T  (K) 

343.15 

Anode  feed  flow  rate,  V F  ( x  10~8  m3  s-1 ) 

2.267 

Geometric  parameters 

Cross-sectional  electrode  area,  As  (xlO-4  m2) 

9 

Volume  of  anode  compartment,  Va  (xlO-6  m3) 

2 

Porosity  of  anode  catalyst  layer,  e 

0.81 

Thickness  of  PEM,  dM  (xlO-6  m) 

200 

Thickness  of  catalyst  layer,  dCL  (xlO-6  m) 

10 

Thickness  of  diffusion  layer,  dDL  (xlO-6  m) 

100 

Design  parameters 

Double  layer  capacity  of  anode,  Ca  ( x  105  F  ur 2 ) 

4 

Double  layer  capacity  of  cathode,  Cc  ( x  105  F  m-2 ) 

8 

Standard  cell  voltage,  U8™  (V) 

1.21 

Standard  pressure,  PSTD  (Pa) 

101,325 

Universal  constants 

Gas  constant,  R  (J  mol-1  K-1 ) 

8.314 

Faraday’s  constant,  F(C  mol-1 ) 

96,485 

Physical  properties 

Molecular  weight  of  water,  Mw  (kg  mol-1 ) 

0.018 

Molar  concentration  of  water,  cH2o  (m°l  m-3) 

55,556 

volume  unit.  Dohle  et  al.  [20]  presented  the  water  content  of  the 
membrane  as  a  function  of  the  relative  humidity,  which  is  defined  as 
the  ratio  of  the  partial  pressure  of  water  to  the  saturated  pressure  of 
water.  Siebke  et  al.  [21  ]  reported  a  more  convenient  equation,  which 
depends  only  on  temperature  and  assumes  that  the  membrane  is 
fully  hydrated: 

if  =  9.38  +  0.138(7  -  273.15)  (16) 


1200 

1000 

800 

600 

400 

200 

0 


(/) 

c 

CD 

Q 


O 


<D 

O 


Cell  Current  Density  [A/m2] 


Fig.  3.  (a)  Dynamic  response  of  cell  voltage  to  step  change  of  cell  current  den¬ 
sity  when  feed  concentration  of  methanol  is  set  to  500  mol  m-3  and  (b)  predicted 
polarization  curves  with  experimental  data  [9]. 
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Table  2 

Summary  of  parameter  estimation  results 


Combining  Eqs.  (15)  and  (16),  the  ionic  conductivity  is  obtained 
as  a  function  of  temperature: 

km  =  (0.07097 -14.88)  exp  (l268  (^  -  1))  (17) 


(a) 


Cell  Current  Density  [A/m2] 
(C) 


Cell  Current  Density  [A/m2] 

(e) 


The  density  of  water  in  Eq.  (10)  depends  on  the  temperature: 

6 

lnpH2o  =  ^aj(T-;rref)i-1  (18) 

1=1 


(b) 


Cell  Current  Density  [A/m2] 

(d) 


Cell  Current  Density  [A/m2] 

(f) 


Cell  Current  Density  [A/m2] 


Fig.  4.  Steady-state  simulation  results  according  to  feed  concentration  of  methanol  from  0.125  to  2  M  with  interval  of  0.125:  (a)  polarization  curves;  (b)  power-current  curves, 
(c)  methanol  crossover  fluxes;  (d)  methanol  crossover  ratios;  (e)  anode  overpotentials;  (f)  cathode  overpotentials. 
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where  the  six  coefficients  are  6.9094,  -2.0146  x  10-5, 

-5.9868  x  10-6,  2.5921  x  10"8,  -9.3244  x  10”11 ,  and 

1.2103  x  10-13,  respectively,  and  the  reference  temperature 
(Tref)  is  273  K. 

The  mathematical  model  involves  six  differential  equations; 
thus  the  same  number  of  initial  conditions  are  required: 

~  i  — init  -  i  — init  ~CL  i  „CL,init 

CMeOHlt=0  -  CMeOH’  cCO2lt=0  -  CC02 »  CMeOF[lt=0  -  CMeOH  ’ 

cS2|t=0  =  c^,  rc|t=0  =  r‘nit,  UceU\t=0  =  U^  (19) 

where  the  superscript  ‘init’  means  the  initial  value.  Since  zero  initial 
hold-up  conditions  are  considered  in  this  research,  the  initial  hold¬ 
ups  of  methanol  and  carbon  dioxide  in  the  anode  compartment  and 
anode  catalyst  layer  will  be  zero. 


in  the  methanol  and  water  tanks,  respectively.  A  sensor  detects  the 
methanol  concentration,  sends  a  signal  to  the  control  valves  con¬ 
nected  to  the  water  and  methanol  tanks,  and  finally  permits  control 
of  the  feed  concentration  from  0  to  2  M.  It  is  assumed  that  the  output 
concentration  of  methanol  (c^ok°Ut)  shows  a  first-order  response 
to  the  input  change  (c^qh11)  *n  a  piecewise  constant  fashion: 


j  Feed, out 

ULMeOH  ,  ^Feed,out  _  ^Feed,in 

df  +  cMeOH  ~~  ALMeOH 


(20) 


where  r  is  the  time  constant  and  I< is  the  steady-state  gain.  The  feed 
pump  supplies  the  cell  with  methanol  solution  at  a  constant  flow 
rate.  The  operating  conditions  and  parameters  are  summarized  in 
Table  1. 


2.4.  Parameter  estimation 


23.  Target  system 

To  describe  the  dynamic  operation  of  a  DMFC,  the  system  that 
consists  not  only  of  the  cell  but  also  the  accessories,  such  as  the  feed 
pump,  storage  tanks  and  control  valves  with  the  methanol  sensor, 
as  shown  in  Fig.  2.  Initially,  2  M  methanol  and  pure  water  are  stored 


Time  [s] 


Fig.  5.  Transient  behaviour  of  DMFC:  (a)  scenario  for  power  density  load;  (b) 
input  and  output  concentrations  of  methanol  feeding;  (c)  concentration  profiles 
of  methanol  and  carbon  dioxide  in  anode  compartment  (C)  and  in  anode  catalyst 
layer  (CL);  (d)  reaction  rates  and  methanol  crossover  flux  (left  axis)  and  methanol 
crossover  ratio  (right  axis);  (e)  cell  voltage  profiles. 


To  verify  the  model,  five  parameters  were  selected  from  the 
model  equations,  namely:  the  charge  transfer  coefficients  of  the 
anode  and  cathode  (aa,  ac),  the  reaction  rate  constants  of  the 
anode  and  cathode  (fca,  kc\  and  the  mass  transfer  coefficient  (/<LS). 
The  parameter  estimation  was  performed  in  the  same  way  as  the 
experimental  method  shown  in  Fig.  3a.  It  minimizes  the  difference 
between  the  experimental  data  and  the  cell  voltage  responses  to 
the  piecewise  constant  inputs  of  the  cell  current  density  to  find  the 
best-solution  set  of  selected  parameters.  The  optimization  model 
for  the  parameter  estimation  can  be  formulated  as 


mm 

71;  e  ft 


I<  J 

EE 

k= 1  j= 1  L 


m(o in 


(Uce\\,jk  U cell  Jk ) 


G 


jk 


s.t.  /(x(t),x(t),y(t),  v)= 0, 


/(x(0),  x(0),  y(0),  v)=0, 

7t\  <  71  [  <  7T2U  (21) 

here  77  (the  set  of  parameters  to  be  estimated)  =  {aa,  oc c,  /<a,  kc,  kls} 
and  all  of  parameters  have  their  own  boundaries;  L  and  U  denote 
the  lower  and  upper  boundaries,  respectively.  The  symbols /and  / 
represent  the  model  equations  discussed  in  the  previous  section, 
and  the  initial  conditions  for  differential  equations.  Each  experi¬ 
ment  k  has  Jk  piecewise  constant  inputs  of  cell  current  density.  ajk 
is  the  variance  of  the  jth  step  input  of  the  cell  current  density  in 
experiment  k.  For  a  variance  model,  the  heteroscedastic  variance 
model  was  chosen,  i.e., 

a2=®2(Uc2ell  +  e)y  (22) 

where  s  is  the  absolute  tolerance  and  has  a  fixed  value  of  1  x  10-5. 
Two  parameters,  go  and  y ,  will  be  determined  from  the  optimization 
under  their  lower  and  higher  boundaries.  To  estimate  the  selected 
five  parameters,  we  used  the  estimation  module  in  gPROMS  [22,23 ] 
and  two  sets  of  experimental  data  for  methanol  feed  concentrations 
of  500  and  2000  mol  m-3  from  Sundmacher  et  al.  [9].  The  polariza¬ 
tion  curves  from  the  estimation  results  show  good  agreement  with 
the  experimental  data,  as  seen  in  Fig.  3b.  The  R2  values  of  the  two 


Table  3 

Six  scenarios  generated  for  the  case  study 
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cases  are  0.946  and  0.902,  respectively.  The  parameter  estimation 
results  are  summarized  in  Table  2. 

3.  Model  verification 

3.1.  Steady-state  results 

Fig.  4  presents  the  steady-state  cell  performances  based  on  the 
feed  concentration  of  methanol.  Six  variables  are  plotted  against 
cell  current  density:  (i)  cell  voltage,  (ii)  cell  power  density,  (iii) 
methanol  crossover  flux,  (iv)  methanol  crossover  ratio,  (v)  anode 
overpotential,  and  (vi)  cathode  overpotential.  Each  figure  thus  has 
16  graphs,  and  they  match  to  the  feed  concentration  of  methanol 
from  0.125  to  2  M,  with  intervals  of  0.125,  increasing  in  the  direction 
of  the  arrow.  These  results  are  obtained  from  the  dynamic  simula¬ 
tions  using  gPROMS  under  potentiostatic  conditions,  by  collecting 
data  when  the  simulations  reach  steady  state.  As  shown  in  Fig.  4, 
the  model  is  well  matched  with  the  steady-state  trends  previously 
reported  in  the  literature. 

3.2.  Dynamic  behaviour 

To  verify  the  dynamic  behaviour  of  the  DMFC  over  the  entire 
operating  procedure,  a  dynamic  simulation  was  performed  from 
start-up  to  shut-down  operation.  The  total  operating  time  was 
1300  s:  the  start-up  (feeding)  for  50  s,  the  normal  operation  for 
750  s,  and  the  shut-down  for  500  s  (long  enough  to  consume  all 
methanol  in  the  cell).  Fig.  5a  and  b  shows  the  strategies  of  the  power 
density  load  and  the  methanol  feeding.  At  start-up  operation, 
no  power  is  loaded,  but  1  M  methanol  is  introduced,  and  nei¬ 
ther  methanol  feeding  nor  power  loading  exists  at  the  shut-down 
period.  During  normal  operation,  300  Wm-2  of  power  density  was 
loaded  constantly,  and  three  different  concentrations  (1.0,  0.5,  and 
1.5  M)  were  fed  sequentially  for  an  equal  period  (250  s)  to  check 
the  effect  of  methanol  feeding  on  cell  performance.  Fig.  5c  shows 
the  concentration  profiles  of  methanol  and  carbon  dioxide  in  both 
the  anode  compartment  and  the  anode  catalyst  layer.  The  methanol 
concentration  tends  to  vary  with  the  feed  concentration;  whereas 
the  concentration  of  carbon  dioxide  shows  a  tendency  similar  to 
the  anode  reaction  rate  shown  in  Fig.  5d.  As  soon  as  the  cell  enters 
shut-down  operation,  the  anode  reaction  rate  drops  dramatically. 
By  contrast,  the  methanol  crossover  flux  decreases  gradually  and 
results  in  a  steady  fall  in  methanol  concentration  in  both  the  anode 


Fig.  6.  Summary  of  54  optimization  results;  six  scenarios  and  nine  optimizations 
per  scenario  varying  NCA. 


compartment  and  the  anode  catalyst  layer.  This  trend  indicates 
that  residual  methanol  remaining  after  normal  operation  results 
only  in  crossover,  which  is  a  waste  of  fuel.  The  anode  reaction  that 
occurs  during  the  normal  operation  varies  slightly  with  the  con¬ 
centration  of  methanol.  On  the  other  hand,  methanol  crossover 
depends  strongly  on  the  concentration  of  methanol  in  the  anode 
catalyst  layer  and  exceeds  the  anode  reaction  when  more  methanol 
is  present  than  required,  especially  when  the  cell  is  in  the  start-up 
or  the  shut-down  periods.  Consequently,  the  methanol  crossover 
ratio,  which  is  defined  as  the  ratio  of  the  methanol  crossover  rate 
to  the  sum  of  the  methanol  crossover  rate  and  the  anode  reac¬ 
tion  rate,  is  a  maximum  in  the  start-up  and  shut-down  periods, 
as  demonstrated  in  Fig.  5d.  The  cell  voltage  drops  suddenly  at  the 
early  stage  of  normal  operation,  and  varies  with  the  feeding  strat¬ 
egy;  it  also  rises  dramatically  to  reach  the  open-circuit  value  in  the 
shut-down  period  (Fig.  5e).  From  the  dynamic  simulation  results, 
it  is  found  that  the  transient  behaviour  of  the  DMFC  depends  on 
the  feed  concentration  of  methanol  and  that  the  performance  of 
the  cell  can  be  improved  considerably  by  reducing  the  amount  of 
methanol  crossover,  that  is,  by  optimizing  the  strategy  of  methanol 
feeding. 

4.  Dynamic  optimization 

4.1.  Model  formulation 

The  aim  of  the  dynamic  optimization  of  this  research  was  to 
find  the  optimum  operating  strategies  that  maximize  the  perfor¬ 
mance  of  the  DMFC  under  a  given  power  density  load.  The  energy 
efficiency  can  be  a  candidate  for  the  objective  function  to  be  max¬ 
imized: 

fnf  ^celD’cell^S  dt 

Energy  =  ^  ^  ^  *  +  J dt)  (23) 

where  tf  means  the  total  operating  time  of  the  DMFC  system.  The 
numerator  of  Eq.  (23)  equals  the  integral  of  the  power  load,  which 
would  be  constant  if  the  scenario  of  power  load  is  fixed.  Hence,  the 
amount  of  methanol  consumed  by  the  reaction  and  crossover  in 
the  denominator  can  be  used  as  an  alternative  to  the  energy  effi¬ 
ciency.  By  adding  the  amount  of  methanol  remaining  after  normal 
operation  is  terminated,  we  formulate  the  objective  function  for 


30  min  1  hr  2  hr  3  hr  6  hr 

Normal  Operation  Time 


Fig.  7.  Variation  of  methanol  consumption  ratio  with  increase  in  normal  operation 
time;  first  bar  is  base  case,  and  total  amounts  of  methanol  consumed  are  indicated 
by  bars. 
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the  dynamic  optimization,  namely: 

r*  f  rtf 

MMeOH  =  /  raAs  d  t+  n|^je0H/ls  dt 

Jo  Jo 

+(^aLcMeOH  ^aCMeOH)lt=tf  (24) 

here  MMeoH  is  called  the  amount  of  methanol  consumed,  which 
consists  of  the  amount  of  methanol  reacted,  the  crossover,  and 
the  remainder.  The  feed  concentration  is  selected  as  a  piecewise 
constant  control  variable  to  be  determined.  Finally,  the  dynamic 
optimization  model  is  formulated  as 

Min  MMeoH  s.t.  /(x(t),  x(t),y(t),  v)=0, 

rFeed 

MeOH 

J(x(0),x(0),y(0),v)=0, 

Ucell  >  UcLel|  @t  =  tf 

^cell  =  ^cell,s  @  t s  <  t  <  ts_ |_i ,  S  =  0,  1 ,  . . . ,  NS 

0  <  <  2000  mol  m  3  @  0  <  t  <  tf 

(25) 

where  /  and  I  are  the  model  equations  and  the  initial  conditions, 
respectively  discussed  in  Section  2;  NS  means  the  number  of  inter¬ 
vals  for  the  power  density  load. 


Fig.  8.  Transient  behaviour  of  NCA  =  1  of  Scenario  1.  (Subtitles  of  each  graph  are 
same  in  Fig.  5.) 


4.2.  Case  study 

The  purpose  of  the  case  study  was  to  determine  the  influences 
of  (i)  the  power  load  and  (ii)  the  number  of  control  actions  on  the 
performance  of  the  DMFC  system  as  shown  in  Fig.  2. 

First,  six  scenarios  were  generated  that  considered  the  real  oper¬ 
ational  phenomena  of  the  system.  All  scenarios  had  four  steps  with 
a  total  operation  time  of  1860  s;  Step  1  was  the  feeding  step  (60  s) 
and  Steps  2-4  were  the  normal  operation  steps  (600  s  for  each  step). 
At  Step  1  of  each  scenario,  no  power  was  loaded  into  the  cell.  From 
Steps  2-4,  three  different  types  of  power  density  were  loaded  into 
the  cell  without  duplication  for  each  scenario:  high  (460  Wm-2), 
normal  (430  W  m-2 ),  and  low  (400  W  m-2 ).  Each  scenario  produced 
215  Whnrr2  during  the  operation.  Table  3  presents  a  summary  of 
the  six  scenarios  of  the  case  study. 

Second,  the  performance  of  DMFC  system  was  also  presumed 
to  depend  on  how  many  control  actions  were  included.  It  can  be 
expected  that  the  increase  in  control  actions  could  enhance  the 
performance  of  the  system.  To  insure  this  result,  nine  dynamic  opti¬ 
mizations  were  performed  for  each  scenario,  varying  the  number  of 
control  actions,  NCA  e  {1, 4,  7, 10, 13, 16, 19,  25,  31 }.  NCA  =  1  means 
that  methanol  is  fed  optimally  in  a  time-invariant  manner.  It  is 
assumed  that  the  minimal  control  interval  is  60  s  and  that  every 


Fig.  9.  Transient  behaviour  of  NCA  =  7  of  Scenario  1.  (Subtitles  of  each  graph  are 
same  in  Fig.  5.) 
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interval  length  is  the  same,  except  for  the  feeding  step.  The  feed¬ 
ing  step  (Step  1)  is  too  short  to  be  divided,  but  other  steps  can  be 
divided  into  up  to  10  sections.  For  example,  NCA  =  25  means  that 
three  steps,  except  for  the  feeding  step,  are  divided  into  eight  inter¬ 
vals;  the  entire  operation  (1860  s)  is  thus  divided  into  a  60-s  interval 
and  24  segments  of  75-s  intervals. 

Fig.  6  summarizes  the  optimaum  values  of  54  dynamic  opti¬ 
mizations  (six  scenarios  and  nine  optimizations  per  each  scenario). 
When  NCA=1,  the  methanol  consumptions  of  each  scenario  are 
between  23.3  x  10-3  and  23.7  x  10-3  mol.  They  drop  suddenly 
according  to  the  last  load  of  power  density  when  NCA  =  4.  After 
that,  MMe0H  decreases  gradually  as  NCA  increases  in  pairs.  As  can 
be  seen  in  Fig.  6,  the  more  that  the  operation  is  finished  with  a 
high  load  of  power  density,  the  more  the  MMe0H  is  located  at  a 
high  position.  If  two  scenarios  finish  the  operation  with  the  same 
load  of  power  density,  the  scenario  that  starts  from  a  lower  load 
of  power  density  consumes  slightly  less  methanol.  The  reason  why 
methanol  consumption  depends  strongly  on  the  final  load  of  power 
density  is  related  to  the  amount  of  residual  methanol,  the  last  term 
of  Eq.  (24).  To  produce  high  power  density,  the  anode  reaction 
requires  an  equivalent  high  concentration  of  methanol.  If  the  opera¬ 
tion  ends  by  producing  a  high  power  density,  a  considerable  amount 
of  methanol  remains,  finally  compromising  the  objective  function. 


This  dependence,  however,  is  increasingly  weakened  as  the  normal 
operation  time  increases.  This  phenomenon  is  explained  in  Fig.  7  in 
the  form  of  a  stacked  bar  graph  and  shows  the  optimaum  ratio  of 
methanol  consumption  as  a  function  of  the  increase  in  the  normal 
operation  time  from  30  min  (equal  to  the  case  study)  to  6  h,  applying 
the  power  load  of  Scenario  2.  The  amount  of  methanol  remaining, 
which  is  mainly  responsible  for  the  dependence,  decreases  from 
20.5  to  2.1%  as  the  normal  operation  time  increases;  whereas  the 
amounts  of  reacted  and  crossover  methanol  increase. 

The  dynamic  behaviour  of  54  optimization  results  can  be  deter¬ 
mined  by  applying  each  optimum  feeding  strategy  to  the  dynamic 
simulation.  To  compare  the  performance  of  complicated  piecewise 
constant  control  with  the  time-invariant  control,  four  optimiza¬ 
tions  were  selected:  NCA=  1  and  7  of  Scenario  1  (low-ending  power 
load  case),  and  NCA=1  and  31  of  Scenario  5  (high-ending  power 
load  case).  Figs.  8  and  9  summarize  the  dynamic  behaviour  of  the 
optimization  results  of  NCA  =  1  and  7  of  Scenario  1,  respectively.  Due 
to  the  piecewise  constant  input  of  feed  concentration,  the  perfor¬ 
mance  of  NCA  =  7  is  much  better  than  that  of  NCA  =  1 ;  the  methanol 
crossover  and  final  concentration  of  methanol  in  the  cell  are  both 
reduced.  Similar  results  are  compared  for  the  high-ending  power 
load  scenario  in  Fig.  10  (NCA  =1  of  Scenario  5)  and  Fig.  11  (NCA  =  31 
of  Scenario  5).  The  complicated  feeding  strategy  (Fig.  lib)  not  only 


Fig.  10.  Transient  behaviour  of  NCA  =  1  of  Scenario  5.  (Subtitles  of  each  graph  are  Fig.  11.  Transient  behaviour  of  NCA  =  31  of  Scenario  5.  (Subtitles  of  each  graph  are 
same  in  Fig.  5.)  same  in  Fig.  5.) 
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Number  of  Control  Actions 

Fig.  12.  Reduction  in  methanol  consumption  of  each  control  strategy  compared  with 
time-invariant  control  (NCA  =  1 ). 

results  in  fluctuating  behaviour  of  the  variables  but  also  exhibits 
good  performance  (i.e.,  low  methanol  crossover  ratio). 

Fig.  12  presents  the  reduction  in  methanol  consumption  of 
each  control  strategy  compared  with  the  time-invariant  control 
(NCA  =  1).  The  more  the  operation  ends  with  a  high  power  load, 
the  more  it  reduces  methanol  consumption.  Increasing  NCA  leads 
to  increasing  methanol  consumption  in  the  cell.  The  piecewise 
constant  control  of  feed  concentration  reduces  the  methanol  con¬ 
sumption  from  18.53  to  22.31%  compared  with  the  time-invariant 
control  under  Scenario  1,  and  6.07-8.69%  under  Scenario  6. 

5.  Conclusions 

The  optimum  operating  strategies  of  DMFCs  over  the  entire 
operation  including  the  start-up  and  shut-down  have  been  stud¬ 
ied.  A  mathematical  model  that  allows  a  description  of  zero  initial 
hold-up  conditions  is  presented  as  a  set  of  DAEs  and  verified 
with  experimental  data  using  the  parameter  estimation  method. 
Steady-state  simulation  results  properly  demonstrate  the  feed  con¬ 
centration  dependencies  of  key  variables  and  dynamic  simulation 
results  rigorously  explain  the  transient  behaviour  of  DMFCs.  Based 
on  this  model,  a  dynamic  optimization  problem  has  been  formu¬ 
lated  as  a  non-linear  programming  (NLP)  problem.  The  methanol 
consumption,  which  consists  of  the  amount  of  methanol  reacted, 
the  crossover,  and  the  remainder  after  operation,  is  selected  as 
the  objective  function  to  be  minimized.  As  a  decision  variable,  the 
methanol  concentration  of  the  anode  feed  stream  is  considered 
in  the  form  of  a  piecewise  constant  control.  In  the  case  study,  six 
scenarios  of  power  density  load  are  generated  to  reflect  real  oper¬ 


ational  phenomena.  By  varying  the  number  of  control  actions,  54 
cases  of  dynamic  optimization  are  performed  and  compared  with 
each  other. 

From  the  case  study,  it  is  found  that  the  performance  of  the 
DMFC  depends  on  the  number  of  control  actions  as  well  as  on  the 
power  density  load.  If  the  overall  power  generation  is  the  same,  the 
more  the  operation  finishes  with  a  high  power  density  load,  the 
less  efficient  the  system  becomes.  Increasing  the  number  of  con¬ 
trol  actions  reduces  methanol  consumption,  although  this  requires 
a  complicated  control  strategy.  It  is  also  demonstrated  that  opti¬ 
mum  feeding  strategies  reduce  the  feed  consumption  from  6.07 
to  22.31%,  depending  on  the  scenario  and  the  number  of  control 
actions.  In  conclusion,  the  optimum  feeding  strategies  obtained 
from  the  dynamic  optimizations  enabled  DMFCs  to  be  operated 
more  efficiently. 
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